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We suggest a procedure for calculating correlation functions of the local 
densities of states (DOS) at the plateau transitions in the Integer Quantum 
Hall effect (IQHE). We argue that their correlation functions are appropri- 
ately described in terms of the SL(2, C)/SU(2) WZNW model (at the usual 
Kac-Moody point and with the level 6 < k < 8). In this model we have iden- 
tified the operators corresponding to the local DOS, and derived the partial 
differential equation determining their correlation functions. The OPEs for 
powers of the local DOS obtained from this equation are in agreement with 
available results. 
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I. INTRODUCTION 

Under the conditions of low temperature and strong perpendicular magnetic field, a 
two-dimensional electron gas exhibits a striking macroscopic manifestation of a quantum 
phenomenon, namely the quantum Hall effect [|1], 0: the Hall conductivity exhibits quan- 
tized plateaus at well denned multiples of e 2 /h (a fundamental constant). In samples where 
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the role of random impurities (disorder) is more important than electron-electron interac- 
tions, the plateaus occur at integer multiples of e 2 /h giving rise to the so-called integer 
quantum Hall effect (IQHE). It is widely believed that in the absence of a magnetic field, 
all wavefunctions for non-interacting, disordered electrons in two dimensions are localized. 
In the presence of a magnetic field however, a delocalized state occurs at the centre of the 
(disorder broadened) Landau level, with energy E c . As one tunes the electron energy E (by 
varying the magnetic field), through the centre of a Landau level, the localization length, 
£, diverges as £ = \E — E c \~ v , where numerical simulations indicate that v ~ 2.3.Q The 
plateaus with differing a xy are separated by these critical points. A theoretical description 
of these points remains one of the most challenging unresolved problems in the theory of 
disordered systems. 

In the field-theoretic approach to the problem of disordered electrons, the presence of a 
static potential preserves the one-body energy, so the frequency becomes simply a parameter 
of the action and one may consider each frequency separately. This reduces the problem of 
disordered electrons in d spatial dimensions to a d- dimensional Euclidean field theory. The 
natural description of the plateau transitions in the IQHE, should therefore be in terms of 
the critical point of some two-dimensional Euclidean field theory, to which one may apply 
the powerful machinery of conformal field theory (CFT) — see for example the books |3|-0. 
Despite persistent efforts over the last fifteen years, the form of the effective action describing 
the critical point is still uncertain. 

The first field-theoretical description of the disordered Landau level was given by Levine, 
Libby and Pruisken || and Pruisken |7] in the form of a non-linear sigma model with a topo- 
logical term. On the basis of this model, Khmelnitskii suggested a two-parameter scaling 
theory || of the IQHE. In this theory, all the (renormalization group) flow lines merge into 
one of a number of fixed points occurring at integer multiples of the Hall conductance (mea- 
sured in units of e 2 /h) corresponding to the existence of plateaus in the IQHE. In addition, 
the flow-diagram contains unstable fixed points, corresponding to the transition states oc- 
curring between plateaus. These flows were justified in a dilute instanton gas approximation 
by Levine, Libby and Pruisken || and independently by Knizhnik and Morozov jlOfl . Un- 
fortunately, the extrapolation towards the conjectured fixed points lies at strong coupling, 
and quantitative results are lacking. In the original derivation of the model, Pruisken and 
collaborators employed the method of replicas, leading to a sigma model defined on the 
manifold U(2n)/U(n) xU(n) with n — > 0. A more rigorous formulation using Efetov's su- 
persymmetry approach [ |TT| , |r2| was given in [I3|, and led to a sigma model on the manifold 



SU(1, 1 |2)/U(1| 1) x U(l|l). The action may be written in the following form, 



S = d 2 x str 



± (^Q) 2 + 1^° (Q [d x Q,d y Q]) + \, 



We note, that in a system of size L, the number of states with wave functions that reach the 
boundaries is ~ p{E)L 2 ~ l / v . Since the density of states, p(E), remains a smooth function of energy, 
this number is always macroscopic. However, the density of delocalized states ~ L~ l / U and goes 
to zero in the limit of infinite sample size. 
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where Q is a 4 x 4 supermatrix satisfying the conditions, 



strQ = 0, Q 2 = 1, (1.2) 

S 3 = diag(l, — 1, 1, — 1) (in the boson-fermion supermatrix representation), a xy (E) is the 
bare Hall conductance at energy E, p is the average density of states at energy E, and r] is 
the imaginary frequency which serves as a symmetry breaking field. At weak disorder, a < 1, 
the bare value of the inverse coupling constant, a -1 , coincides with the bare longitudinal 
conductance a xx (E) at energy E. To be more precise, one can directly relate the bare 
parameters of the action to experimentally observable quantities of a mesoscopic system of 
size L ~ I (the mean free path). At E — E c , one has <r° = |, and the model is expected to 
be critical at large distances. 

The second term appearing in (|1.1|) is topological, and despite the fact that its presence is 
crucial for the critical behaviour, its effect cannot be spotted in a perturbative expansion in 
powers of a; it does not contribute to the equations of motion and hence does not contribute 
to the loop expansion of the beta function. The effects of the topological term become visible 
only for samples of size greater than £ ~ I exp^cr!^ 2 ], where I is the electron mean free path. 
In the model ( |1.1| ), with a xy = 0, the length scale £ corresponds to the localization length; 
with o~ xy = 1/2, this scale is the transmutation length (in field-theoretic jargon) and signifies 
a crossover to the regime of universal critical fluctuations. 

It is widely believed that the model ( |1 . 1| ) has a non-trivial infrared fixed point at some 
a* and o~ xy = 1/2. However, one cannot simply substitute a* with the experimental value 
of {o~ xx )~ l (the inverse longitudinal conductance) in order to obtain the effective action at 
the critical point. The reason for this is that the fixed point occurs at strong coupling, 
where the fundamental fields of the model ( |1.1| ) are strongly fluctuating. In this regime, 
the coupling constant a, being dependent on the regularization procedure, can no longer be 
identified with any measurable quantity. This phenomenon is well known in asymptotically 
free gauge theories — the gauge coupling constant has meaning only at short distances, and 
cannot be defined in a universal way in the infrared.f] In QCD, for example, gluons do not 
exist as asymptotic states in the infrared — there is a mass transmutation phenomenon 
and the emerging degrees of freedom are massive. In the quantum Hall effect, however, we 
expect the emerging infrared theory to be a massless CFT, and so one should not stretch 
the analogy too far. A better analogy is the Seiberg duality |L3| appearing in N = 1 
supersymmetric gauge theories, with N c colours and Nf flavours, in the so-called conformal 
window, 3N c /2 < Nf < 3N C . There is an infrared (IR) fixed point in this theory, in the 
vicinity of which, the theory can be described as an infrared limit of another gauge theory 
(with the same number of flavours and Nf — N c colours). If Nf is close to 3N c /2 the original 
theory is strongly coupled near the IR fixed point, whereas the second is only weakly coupled; 
the appropriate IR behaviour is given by the second (dual) formulation. 



2 Notable exceptions to the rule are theories with spontaneously broken gauge symmetry (in which 
case zero-charge U(l) can also emerge) where the coupling constant can be defined in the infrared 
through either an effective Fermi four-fermion interaction or Thompson scattering. 
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Another example of this nature occurs in the deformation of the minimal model M } 



v 



of two-dimensional CFT, by the operator $ 13 |15|,|16fl. The theory in the infrared can be 
described either as a strongly fluctuating asymptotically free M p model, or as a weakly 
fluctuating infrared free M p _i model, deformed by the irrelevant operator $3^. Again, it is 
the second (dual) description which is required in the infrared. 

These analogies encourage us to believe that the theoretical description of the plateau 
transitions in the IQHE will ultimately lie in a reformulation of the model in terms of 

degrees of freedom more appropriate for the infrared region. 

An alternative approach to the study of the plateau transitions is based on the Chalker- 
Coddington Network Model [|T7]]. This model has been reformulated in Hamiltonian form as 



a Replica Spin Chain [0,^9], and as a Superspin Chain |20|-p3|. The Superspin Chain, is 



a lattice model with nearest neighbour antiferromagnetic exchange between 'spins', which 
are generators of the gl(2|2) algebra. The advantage of such reformulations is that one 
may approach the critical point more accurately by tuning the parameters of the network. 
This makes the model indispensable for numerical simulations. However, the Chalker- 
Coddington model (or the spin chain models) are lattice theories, and one still needs to 
derive their continuum limits in order to describe the critical fluctuations. The derivation 
of such limits remains an open problem. 

Exasperated by long and futile efforts to find a rigorous derivation of the effective action 
at the critical point, Zirnbauer has drawn upon a great many sources and conjectured its 
possible form [|4j]PL Such a conjecture ought to be justified, amongst other considerations, 



by its ability to reproduce existing results for physical quantities at the critical point. In p4 
this has been successfully done for the two-point conductance. In the present paper we shall 
concentrate on the correlation functions of the local density of states (DOS). 

The structure of our paper is as follows: in section |TJ we discuss a possible form for 
the correlation functions at the transitions. In particular, we discuss the operator product 
expansion (OPE) for moments of the local density of states p q , which reflects the multifrac- 
tality of the critical wavef unctions. The two-point correlation function of the local DOS is 
given as a particular limit of a four-point function. The latter function includes two opera- 
tors with anomalous dimension zero (vacuum insertions) which serve to redefine the ground 



state of the theory. In section III we discuss more fully the important role played by these 



operators. In section [IV] we discuss a possible connection between the properties of the 



required two- dimensional conformal field theory, and non-perturbative results |27[ obtained 
by studying the one- dimensional limit of model ( |1 . 1| ) . The analysis of the one-dimensional 
case brings us to the idea that one can study the behaviour of the local DOS in terms 
of the SL(2,C)/SU(2) Wess-Zumino-Novikov-Witten (WZNW) model. In sections [V] and 
f\/T\ we study one of the candidate theories for describing the plateau transitions, namely 



the PSL(2|2) WZNW model |24j| . The non-compact bosonic sector of this theory is the 
SL(2,C)/SU(2) WZNW model. In section [VTT| we use the SL(2, C)/SU(2) WZNW model 
to study the correlation functions of the local DOS. Finally we present concluding remarks 



3 Here we would like to mention the recent paper [25] containing interesting suggestions with 
respect to derivation of the critical theory. 
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(including open problems) and technical appendices. 



II. CONSTRAINTS ON CANDIDATE THEORIES 

Any conformal field theory purporting to describe the plateau transitions must satisfy a 
number of constraints, some of which have been discussed by Zirnbauer [ 24(] . For the purpose 
of this paper, we believe the following are important requirements: 

1. The model should be defined over a manifold associated with some supergroup, and 
must contain a non-trivial (non-unit) operator with zero scaling dimension (see, for 
example |TT||). 



2. The model must reproduce all known scaling dimensions. 

3. Restricting our two-dimensional CFT to a finite-width strip geometry (in a suitable 
limit to be discussed below, and in section |TV|) one should recover certain known 
quasi-one-dimensional results. 

We shall now take each of these constraints in turn, and consider some of their consequences. 



A. Supermanifold and Zero-Dimension Operator 

Given our belief that a (presently unknown) non-perturbative treatment of model ( |1 . 1|) 
will ultimately describe the plateau transitions in the IQHE, it is entirely natural to expect 
the resulting effective action to be defined over some supermanifold. However, in order 
to combine conformal symmetry with bosonic and fermionic degrees of freedom, one will 
ultimately face a zero-mode problem: if the effective action at the critical point contains 
only terms with derivatives, the integral over the zero modes (constant fields) leaves the 
partition function ill-defined. Explicitly, 



dQ = J dB J dF (2.1) 

is indeterminate since the integral over the fermionic fields, F, is zero, whereas the integral 
over the bosonic fields, B, is infinite (assuming that manifold is non-compact as in the 
original model ( |1 . 1| ) ) . In the sigma model ( |1.1[ ) this problem is solved by the presence of the 
additional term 

S v = rj J d 2 xstr(£ 3 Q) (2.2) 

This term is responsible for energy level broadening in the original theory. The correspond- 
ing operator has zero scaling dimension^ and its value does not change in all orders of 
perturbation theory. The term is highly relevant and generates a length scale L v ~ rj~ l l 2 . 



4 We note that the present theory is a purely spatial one, and all dimensions are expressed in 
terms of units of length, not time. From this point of view time is dimensionless and the frequency, 
w ~ l/[Time], corresponds to a dimensionless operator. 
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The zero modes are also eliminated when one considers an open system. The possibility of 
escape through the boundaries gives rise to non-uniform (i.e. energy-dependent) broadening 
of the energy levels. For this reason closed and open systems may require slightly different 
approaches. In this paper we consider only closed systems and shall not discuss any problems 
related to the conductance. In a closed system one has to introduce boundary conditions as 
was suggested by Zirnbauer |26| and subsequently used by various authors |27|,|28|. In the 



field theory approach these boundary conditions can be realized by the insertion of certain 
operators into correlation functions. Such insertions naturally lead to a dependence of the 
correlation functions on the distance Lb from the boundary. However, if the boundary 
operator contains a non-trivial zero dimension operator (we shall call this operator , one 
can obtain a finite answer for the correlation function on sending the boundaries to infinity. 



B. Local DOS 

In a closed system one cannot study conductance. The only remaining quantities to study 
are the wave functions. Since wave functions are not observables, the original sigma model 
does not generate them directly, but rather allows one to calculate correlation functions 
of the local density of states (DOS). We believe that these correlation functions can be 
expressed solely in terms of the (non-compact) bosonic sector of the full theory. (For the 
one-dimensional case this point is discussed in [O], and for two dimensions in ||28|| .) 

In a system of finite size where all energy levels are discrete, the local DOS is defined as 

p(x,E) = J2\Mx)\ 2 S(E-E a ) (2.3) 

a 

where a denotes eigenstates of the system. However, this definition is not suitable for 
calculation of averages of the local DOS because products of delta functions yield infinities. 
These infinities are removed if one introduces a finite level broadening replacing the delta 
functions by the Lorentzians: 

5{E-E a )^- JL 2 (2.4) 

7T (E - E a y + if 

As we have mentioned above, in the sigma model the level broadening comes from term 
Q2.2| ). This term does not break conformal invariance provided the level broadening is much 
smaller than the mean level spacing, rjL 2 << 1. Even such small broadening does the trick of 
removing infinities from correlators of the local DOS. In what follows we consider normalized 
powers of the DOS defined as 

W = (VP) q (2-5) 

In the context of the sigma model ( p..l|) the normalized [p] q is identified with {rjsti^Qk)} 11 
where k is a certain constant matrix. 

In a closed system the wave function on the boundary does not depend on the fields, that 
is equal to unity. However, it may occur (and we shall further explain why this is the case 
for the theories in question) that the unit operator is not a primary field of the problem. 
Then one needs to decompose the unit operator in terms of the primaries: 
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I(z) = + / dfi(p)%(z) 



(2.6) 



In this decomposition nas zero conformal dimension and the other fields ^ p (z) have 
positive dimensions 0. 

Working with the normalized DOS and the boundary operators present one can define 
the correlation functions in the infinite system (we imagine this system as a long strip of 
length 2L with periodic boundary conditions in the transverse direction) ). In the limit 

lim (I(-L)( V py(x)( VP y(y)I(L)) = lim (* (-L)W(x)\ffi\(y)ML)) (2-7) 

V L 2 ^0 

only \l/o operators survive as L — > oo. Taking this into account and also that expression 
( ^.7[ ) is valid in the strip geometry, one obtains the following expression for the two-point 
correlation function in the infinite plane: 

F»1M= , lim, (\fl(Ti)\fl(r a )M**)M**)) ( 2 - 8 ) 

|ra|— >0,|r4|— >oo 

Here we use the overbar to denote disorder averaging, and use the angular brackets to 
denote the correlation functions of our field theory. In the case of strip geometry which can 
be obtained from the plane by a conformal transformation, the point of origin is mapped 
onto minus infinity. We shall return to a more detailed discussion of this procedure towards 
the end of Section IV. 



C. Scaling Dimensions 

The theory we seek should reproduce all known scaling dimensions associated with the 
plateau transition in the IQHE. The most famous of these is the localization length exponent, 
v, mentioned in the introduction.^ 



5 An analogue of this decomposition is equation (3.7) of 



W W {X,t) = 2X 1 ! 2 



Kl (2X^) + - r^ T ^sinh^^ p (2X 1 /2 )exp (_i±Z! r) 



where K n is the modified Bessel function of the third kind and (X, r) is a solution of a 
differential equation 

8 2 „\ T „ W , 8WW(X,t) 



dr 



with the boundary condition W^(X, 0) = 1; here r plays the role of z. 
We see that at zero r this is indeed the decomposition of unity and at infinite r only the first 
term 2X 1 I 2 K 1 (2X 1 / 2 ) survives - which is an analogue of ^o(z). 

6 If we assume that the correlation length is generated by a single operator, of scaling dimension d 
coupled to (a xy — 1/2) and assume further a linear relationship between the Hall conductance and 
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Other available information comes from the behaviour of the local DOS, p(E,r), in a 
two-dimensional sample of side L. At the critical point, p(r) = p(E = E c , r), they are known 
to satisfy the following fusion rulesQ | f2T^ : 

p p {v l )p q {v 2 ) ~ |n - r2 \-d(p)-d(q)+2d(p+ q ) pp+9( r2 ) (2.9) 

p^)~L- d{q) (2.10) 

where, once again, the bar stands for disorder averaging; we reserve angle brackets for field- 
theoretic correlation functions. (Recall that the local DOS are normalized in such a way that 
their correlation functions are finite in the limit of vanishing level broadening.) Numerical 
simulations (see POI] and references therein) indicate that[| (see Appendix |A 3| ) , 

d{q)=2q(l-q)/k, 2/k = 0.28 ± 0.03. (2.11) 

This gives 6 < k < 8. 

The quadratic dependence of d(q) gradually becomes linear for \q\ > 2.5. We note that 
the scaling dimensions become negative for q < and q > 1, well inside the interval of validity 
( |2.11| ). This suggests that the conformal field theory we are looking for is non-unitary. 



D. Quasi-One-Dimensional Results 

Transfer matrix calculations for model ( |Llp in a quasi-one-dimensional, closed (no ex- 
ternal leads), infinite sample, yield the following form for the two-point correlation function 
of the local DOS (see equations (3.58), (3.63) and (3.64) of |27|): 

RNMW = / rfMp)l(0|Q 9 b)rexp(-|a; 12 |d(p)/4^, (2.12) 



where 



d(p) = l+p 2 (2.13) 
, , 2psmh.(irp) 

dp{p) = — f-^dp 2.14 

poo 

(0\Q"\P)= / dX X q - 2 W (X)W^ 1+ip)/2 (X) (2.15) 
Jo 

^\T[ q +(l + zp)/2]\ 2 \T[ q -(l + zp)/2}\ 2 



\E — E c \ in the vicinity of the transition, we conclude that the scaling dimension of the operator is 
d = 2- l/v~ 1.57. 

7 We adopt the usual CFT nomenclature in anticipation of the Virasoro-primary nature we shall 
subsequently motivate for the local DOS. 

8 The proportionality constant in this expression is denoted as 2/k with some hindsight: the 
parameter k will ultimately be identified as the level of the SL(2,C)/SU(2) WZNW theory. 
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and £ is the correlation length. The functions Wj satisfy the following eigenvalue equation 

~ X ) Wl{X) = 3(3 + ^ W iW- (2 - 16) 

The solutions to this equation may be written in the form 

Wj{X) = 2VXK {1+2j) (2Vx) (2.17) 

where K p is the modified Bessel function of the third kind. The relevance of these quasi- 
one-dimensional calculations for the critical theory in two dimensions, will be discussed in 
detail in section |TV| . We simply note here that the correspondence emerges when one places 
the CFT on a strip of width 2nR. If the infinite plane is parametrized by the coordinates 
(z,z), and the strip by (w,w), one may effect this mapping by means of the conformal 
transformation w = Rlnz. Under certain conditions (in particular, one has to consider 
R ~ £) such a strip may be regarded as a quasi-one- dimensional wire. 



III. DISCRETE AND CONTINUOUS REPRESENTATIONS 



To the reader familiar with conformal field theory, equations (|2.9| ) and (|2.12| ) may appear 



contradictory. If p q is a primary field, under what representation of the symmetry group of 
the theory does it transform? Whatever the symmetry group of the 'grand' theory turns 
out to be, we expect it to be some subgroup (or coset) of Gl(2|2), and as such may support 
both discrete and continuous representations. Equation fl2.9|) suggests that p g transforms 
according to a discrete representation, whereas equation ( |2.12| ) suggests that p q is not a 



single primary field, but a linear combination of fields with different scaling dimensions.^ 

We suggest a resolution of this paradox which is independent of any particular form of 
the plateau theory. The solution is intimately connected to the zero-mode problem and the 
existence of a (non-trivial) zero- dimension operator discussed in section [II A . We suggest 



that the two-point function of the normalized (see the discussion around Eq. ( |2.5p ) local DOS 
is given by the limit of the four-point function of some CFT, 



M(n)M(r 2 )= Km ([^(rO^]^)^^)^^)). (3.1) 

|r 3 |— ►0,|r 4 |— >oo 

Let us discuss this suggestion in detail. Invariance under the projective transformations 
of the plane[|], restricts the above four-point function (which we denote by G q (l, 2, 3, 4)) to 
have the form 

G < ,(l,2,3 ) 4) = r ^ 5r ^,(z,z); z = (3.2) 

212 4 9 -231-242 



9 A simple resolution of the paradox that the integral (2.12) generates simple powers is denied by 
the analytical properties of the exponential function which prevent one from deforming the contour 
of integration onto the poles of the matrix element ( |2.15| ). 



10 Namely transformations of the form w = (az + b)/(cz + d), where ad — be = 1. 
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where we have used the fact that p q has conformal dimension h q {= d q /2), and \&o has zero 
conformal dimension. By definition, the function ( |3.2|) is invariant under the permutation 
of Z\ and z 2 which implies 

F q (z,z)=F q (l/z,l/z). (3.3) 



In order to reproduce the fusion rules of the local DOS ( |2.9| ) one needs to consider the limit, 
z\2 ~~ > 0. The three-point correlation functions are fixed by conformal invariance, their 
holomorphic dependence being given by 



(O hl (l)O h2 (2)O h3 (3)) = C 12 3 z^~ h ^z- x t- h3+h2 z^- hM (3-4) 

where C123 are the so-called structure constants of the theory. Using this fact, one obtains 
the desired limit 



(^](l)[^](2)* (3)*o(4)> - \z 12 \-* h « +2h " ([p 2 «](2)M/ (3)v|/ (4)> 

-fOO 
'2q 



Zl2 \-±h q +2h 2q C W J ^ J -2h 2q ^ ( 35 ) 



(recall that z 3 = 0, z 4 — ► 00). This result fixes the asymptotics of T q {z) as 1 — z — > in the 
following manner: 

F q (z) ~ (l-z) h *«. (3.6) 

As we shall now demonstrate, the two-point function of the local DOS in the strip geometry 
explores different asymptotics of T q {z). One can map the plane to the strip by the trans- 
formation w = Rlnz; point 3 goes to —00 and point 4 goes to +00. As a result one obtains 
(the holomorphic part of) the correlation function as follows: 

G 9 (l,2,3,4) = [2Rsmh(w 12 /2R)]- 2h «F q [exp(w 12 /R)} . (3.7) 

Notice that in the strip geometry the correlation function is translationally invariant. In the 
limit Re w 12 ^> R the behaviour of this function is governed by the asymptotics of J r q {z) for 
\z\ ^> 1. In the limit Re^2i 3> -R it is governed by the asymptotics for \z\ 1. (These two 
limits are related by the crossing invariance condition (|3.3| ).) We see that the fusion rule 
(|2.9|) and the expansion in continuous representations (|2.12|) appear in different channels of 
the four-point function. 

Thus the two-point function of the local DOS explores different asymptotics of the func- 
tion T in different geometries. In this way one may to expect to resolve the paradox. 

We conclude this Section with a list of Operator Product Expansions (OPE) which are 
necessary to reproduce the above results. First we shall specify our normalization conven- 
tions. We normalize all two-point correlation functions of primary fields Oh from discrete 
representations as follows: 

(O h {l)O hl (2)) = 6 hth ,zg h . (3.8) 
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The operators from continuous representations V p (z) are normalized by the invariant mea- 
sure dp, /dp on the groupQ: 

^ v -^ = (3 ' 9) 

Adopting these conventions, one may write the OPEs in the following manner: 

[p 9 ](l)[p p ](2) = C { q q p +P) \z 12 \~ 2h "- 2h ^ +2h ^) [p^}(2) + • • • (3.10) 

[/](l)*o(2) = \z 12 \-^C^ (2) + J d^p)\z 12 \-^ 2h ^Cl q {p)V p {2) + ... (3.11) 

where the ellipsis stand for less singular terms. With these conventions one finds the following 
asymptotics for T{z): 

F(z^0) = C%C 2 «\z\ 2h ^ + --- (3.12) 



T{z — > oo) = \z 



[CU 2 + J dp(p) \C$J 2 \z\~ 2h ^ + ■••]. (3.13) 



IV. CFT ON THE STRIP 

A. General Considerations 

In the absence of a rigorous derivation of the critical model describing the plateau tran- 
sitions in the IQHE, one is forced to make some assumptions about the general form of the 
action. Given the form of the model ( |1.1| ), it is natural to assume that this action is of the 
sigma model type, and probably includes the Wess-Zumino term. The reason for including 
the Wess-Zumino term is related to the fact that nearly all critical sigma models that we 
know of require such a term to ensure criticality. Models of this kind are known as Wess- 
Zumino-Novikov-Witten (WZNW) models; their actions either have a full group symmetry, 
G, or are defined on a coset space, G/H, with H being a subgroup of G. The symmetry 
manifold of the required WZNW model is some supergroup manifold, whose symmetry is 
almost certainly greater than that of the original model ( |1 . 1[ ) . 

The WZNW action may be written in the form 



S = j d 2 x [^/gg^G ab [X)d fl X a d u X b + e^B ab [X}d fi X a d u X b ] 



(4.1) 



where X a are fields representing the coordinates on some group (or coset) manifold, G ab is 
the metric tensor on this manifold, and B ab is an antisymmetric tensor. We define the action 
on a curved (world sheet) surface with a metric g^ u (fi, v = 1, 2). 



n It is highly likely that for the true theory, the measure coincides with the one given by equation 
( 2,14| ) and h(p) = c(l +p 2 )/8. However, in this general discussion we do not need to be so specific. 
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A very important feature of the action (|4.1|) is that the second term does not contain 
the world sheet metric. Consequently, the classical stress-energy tensor, T^ u = SS/Sg^, is 
determined solely by the first term. In particular, the most important components for the 
critical model are given by, 

T zz = G ab [X]d z X a d z X\ T- z - z = G ab [X}d- z X a d- z X b (4.2) 

where z = xq + ixi, z = xo — ix\. Here, the reader should not get the false impression that 
the Wess-Zumino term is not important. Unlike the topological term in action ( |1.1|) it does 
contribute to the equations of motion. Since the model ( |4.1| ) is supposed to be critical, these 
equations (to be understood as identities for correlation functions in the quantum theory) 
are 



T zz = 0, d- z T zz = 0, d z T- z - z = 0. (4.3) 

Their fulfillment depends on the Wess-Zumino term through the dynamics of the underlying 
fields X a . 

The smallness of 1/k means that there are many fields in the theory with conformal 
dimensions much smaller than unity. This gives weight to the idea that the critical point 
occurs in the region where the coupling constant of the sigma model is relatively small. Thus, 
one may attempt to describe the critical point using the semiclassical approximation. We 
shall formulate the semiclassical approximation with the specific aim of establishing contact 
between our calculations, and the calculations of correlation functions of the local DOS for 
quasi-one-dimensional systems performed by Mirlin |27p and Fyodorov and Mirlin |[32|| . 

In the infinite plane, parametrized by coordinates (z,z), conformal invariance restricts 
the two-point function of primary fields of conformal dimension (h,h) to be of the form 

*i)0(z2, z*)) = (zi - z 2 )- 2h ( Zl - z 2 )- 2h . (4.4) 

Under a conformal transformation of the plane, w = w(z), this correlation function trans- 
forms like a tensor of rank (h, h): 

2 / dz \ ^ / dz \ 

(0(tWi,tDi)0(w 2 ,w 2 )) = JJf— J (-7^) zi)<t>{z 2 , z 2 )). (4.5) 

i=i V / Wi \ / Wi 

One may pass from the infinite plane to a strip of width 2ttR by means of the conformal 
transformation w = Rlnz. Combining this transformation with (|4.4|) and ( |4.5|) , one obtains 
the two-point function in the strip geometry: 

((l)(w l ,iD 1 )(f)(w 2 ,iD 2 )) = [2Rsinh{w 12 /2R)]- 2h [2Rsmh{w 12 /2R)]~ 2h . (4.6) 

Introducing coordinates (r,cr) along and across the strip respectively (w = r+ia, w = r—ia; 
— oo <r<oo, < a < 2nR ), one may expand (|4.6|) in the following manner, 



( ( p(r 1 ,a 1 ) ( f ) (T 2 ,a 2 ))=J2 J2 C ™ e^ h+h+n ^/ R e -^-h+m)M/R (4 . 7) 



n=0 m= 
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One may also obtain the two-point function in the operator formalism, leading to the 
Lehmann expansion: 



(0(r 1 ,a 1 )0(r 2 , ( T 2 )) = ^|(O|0|a)| 2 e-^^l-^l- 2 l (4.8) 

a 

where E a and P a are the eigenvalues of the Hamiltonian and the momentum operator 
respectively, in the state \ot). Comparing ( |4.7| ) and ( |4.8| ) one obtains a relationship between 
the eigenvalues of the Hamiltonian and the momentum operator, and the scaling dimensions 
in the corresponding CFT: 

„ h+h+n h—h+m 

E a = , P a = . (4.9) 

Restricting our attention to fields with h — h — d/2, one may rewrite ( f4.6| ) in the form 

Mn, <7iMr 2 , <r 2 )) = R- 2d [2 cosh (t 12 /R) - 2cos(a 12 / R)}~ d . (4.10) 

One observes that for r 12 ^> R the asymptotic form of the correlation function is independent 
of a, 

(0(ti,<7i)<Kt 2 ,<7 2 )> ~ R~ 2d exp {-dr 12 /R) {m > R) (4.11) 



and in this limit one should set n = m = in equations (|4.7[) and ([4.9|) . 

Let us now place the model ( |4.1|) on a thin strip of width 2ttR, and neglect any a- 
dependence of the fields X a . From our considerations in the previous paragraph, such a 
procedure preserves the (large r) asymptotics of the correlation functions. The action ( |4.1|) 
becomes 

S = 2nR J drG ab [X]d T X a d T X b , (4.12) 

which may be recognised as the action for a free, non-relativistic particle, of mass m = AttR. 
The corresponding Hamiltonian is the Laplace-Beltrami operator (multiplied by — l/2m), 

k -^hM^^) (413) 

As we have already established, the eigenvalues of this Hamiltonian are related to the 
spectrum of scaling dimensions in our CFT by equation (O). Moreover, solution of this 



Schrodinger equation allows one to obtain explicit expressions for the eigenstates, \a), ap- 
pearing in 



B. Emergence of the SL(2, C)/SU(2) symmetry 



We now wish to establish contact with the one-dimensional calculations based on the 
original model ( |1.1| ). As is well established, the scaling trajectories for different values of 
cr° only start to deviate at the scale £ ^> / (the mean free path). For the sigma model this 
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scale is in the deep infrared, whereas for the (unknown) critical theory this scale serves as an 
ultraviolet cut-off. The critical theory is conformally invariant, and one may map it to the 
strip in the manner described above; we assume that the general form of such expressions 



holds for strips as narrow as R ~ £. Similarly, we assume that the results |^7| hold for strips 
as wide as £. 



On the one hand, the functions ( 2.17Q entering into the matrix elements ( 2.15 ) should 



be eigenfunctions of the Laplace-Beltrami operator for the critical model we seek. On the 
other hand, as we are going to show, these functions are solutions of the eigenvalue problem 
for the Laplace-Beltrami operator on the manifold SL(2, C)/SU(2). This is striking because 
the SL(2,C) symmetry was not the classical symmetry of the original sigma model and in 
fact appears due to the elaborate limit 77 — described in the previous sections. 

Let us describe the details. An arbitrary element h G SL(2, C)/SU(2) admits the 
following decomposition (see equation (18) of [31] and make the trivial replacements 
(0, //+,//_,) -> (6>,/i, /!*)), 

g c; c ?) ^ 

where 9 6 R, n G C, and fi* is the complex conjugate of \i. Adopting this parameterization, 
one may write the SL(2, C)/SU(2) WZNW action in the following form (see equation (19) 

of 101), 

S = A J d 2 x [4d6d6 + e 2e dfxdfx*] . (4.15) 

One may now read off the corresponding metric G a b (c.f. ( (4.1|) ) and deduce the form of the 
Laplace-Beltrami operator (|4.13| ) on the SL(2, C)/SU(2) symmetric space. The resulting 
eigenvalue equation reads (X = exp(— 6)): 

~ { X2 ^ + X WdJl) F ^,lS,X) = *F x 0hvf,X). (4.16) 
The functions Wj satisfying equation ( |2.16| ) are related to the following eigenfunctions of 

Fj = exp[ik/2* + ik*/j]Wj(X/\k\) (4.17) 



with eigenvalues A = + 1). At the same time the /i-independent solutions of Eq.( (4.16| ) 



represented by the functions X~ J have the same eigenvalues. We suggest that these functions 
with j = —q represent [p] q operators. The matrix element fl2.15|) is then just a Clebsh-Gordan 
coefficient. 

These facts establish a relationship between the SL(2, C)/SU(2) WZNW theory and the 
results of the one-dimensional calculations described in section II. The functions Wj appear- 
ing in the matrix elements ( |2.15 ) are related to the eigenfunctions Fj of the Hamiltonian 



( |4.16| ). Both the vacuum state j = and the excited states (the states with j = —l/2+ip/2) 
belong to the coherent state representations of the SL(2, R) group (see Appendix [A4j and |33| 
for details). The representations of the excited states are distinct in the respect that their 
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angular momentum is a complex number with continuously varying imaginary part. The 
scaling dimension of the vacuum state is zero whereas those of the excited states are pro- 
portional to A = (1 +p 2 )/4. They coincide with the exponents in equation ( |2.13| ). The 
constant wave function, though being a solution of equation (|4.16| ), is not orthogonal to the 
eigenfunctions and therefore does not belong to the basis of eigenstates. In the subsequent 
sections we shall study the SL(2, C)/SU(2) WZNW model in more detail, and show how it 
may emerge as an independent subsector of some 'grand' (supergroup manifold) theory. 

The appearance of the SL(2, C)/SU(2) symmetry absent in the original sigma model 
where it is explicitly broken at finite frequencies is a highly non-trivial fact and deserves 
comment. The calculations of the local DOS performed in |^7|,|3^] employ the general reg- 
ularization procedure described in Section II. Taking the limit 77 — > and keeping (rjQ) 
constant involves working in the limit of very large values of fields. In this limit the 77-term 
Q2.2|) in the action becomes of order of unity (such that 77-dependence disappears) and con- 
tributes the term linear in X to the Schrodinger equation ( |2.16| ). One can say that this term 
is some sort of "quantum anomaly" . 

Thus, in the original sigma model the term linear in X is generated as a potential 
contribution. On the other hand, as we have seen from equations (|4.16|) and ( 4.17|) , equation 
( |2.16| ) can be interpreted as a subsector of the pure Laplace-Beltrami-type equation ( |4.16[ ) , 
which contains no potential terms. This subsector specifies representations of the particular 
type (|4.17| ). It is important that the representations of [p] q (that is the discrete series) 
are solutions of the same equation ( |4.16|) , but in the //-independent subsector. This fact 
allows one to consider (^4.17|) and [p] q as representations of the same group. Thus the one- 
dimensional limit has a hidden dynamical symmetry and this fact gives another argument 
in favour of our conjecture that the critical point possesses the SL(2, C)/SU(2) symmetry. 



V. GAUGED WZNW MODELS 



In this section we consider the WZNW model on the supergroup manifold 



PSL(2|2)=SL(2|2)/GL(1) (see Appendix |X| for more details on supergroups). These two- 
parameter models are quite remarkable in that they are believed to be conformal along 
(marginal) lines of parameter space, rather than isolated points |34],|35|. Zirnbauer has pro- 



posed a theory describing the plateau transitions in a region of this parameter space not 



endowed with the Kac-Moody symmetry [24]. For the purposes of this paper we consider 
the model where the Kac-Moody symmetry is present (the 'Kac-Moody' point). In par- 



ticular, in section [VI], we demonstrate that the bosonic and fermionic sectors decouple; the 
non-compact bosonic sector is described by the SL(2, C)/SU(2) WZNW model. 

The standard Goddard-Kent-Olive (GKO) procedure for dealing with such coset 
spaces may be formulated as a gauged WZNW model (for more details see 0-§). The 
equations for correlation functions can be obtained by the gauge dressing of the conventional 
Knizhnik-Zamolodchikov equations J37j. We consider the WZNW model on the SL(2|2) 
group first and then gauge away the GL(1) subsector. In order to make our arguments 
transparent, we discuss some general aspects of the theory of WZNW models. 

A convenient starting point for discussing the coset G/H, is the (Euclidean) WZNW 
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action on the supergroup manifold G, 



S[g] = S [g] + kT[g}, (5.1) 

where 

So[ 9 ] = ^ J ft-fteiPg^dtf) (5.2) 

t[g] = d 3 xe^str {g-%g g-%9 g-%9) . (5.3) 

The matrix field g(£) is taken to be an element of some supergroup G, the hats indicate 
the presence of the supertraceQ, £ M = (£\£ 2 ) are the coordinates of our two-dimensional 
(Euclidean) space, A 2 and k are dimensionless parameters (H — 1). The Wess-Zumino term, 
T[g], is defined by the integral over the three-dimensional ball with coordinates x a ; the 
boundary being identified with our two-dimensional space. Since the Wess-Zumino term 
contributes to the equations of motion (see below), and hence to the perturbative beta 
function, it is not to be confused with the topological term. 

The action (|5.1| ) satisfies the Polyakov-Wiegmann identity^ f38|. 

d 2 £^ u str {a- 1 d^abdvb- 1 ); = — - ——. (5.4) 

2a o7T 

The classical equation of motion follows from the requirement that the action be stationary 
(0 = 5S = S[g + Sg] — S[g]) under the replacement g — > g + 5g = g(l + g^Sg). Substituting 
the latter form of the variation into (|5.4|), and keeping terms of order 5g, 



= 5S = -J dH^ u M9- l d,gd v {g- l 5g)) d v {u^g-%g) = 0. (5.5) 

Thus for A 2 = 4n/k the field equation becomes 

dig^dg) = ^ didgg- 1 ) = 0. (5.6) 
We shall denote the action S[g] at the Kac-Moody point (A 2 = 4ir/k) by the symbol VP[g]: 

W[g] = kf[g] = ±-J d^st^g-^g) + kf[g]. (5.7) 



12 In the subsequent analysis removal of the hats in any term corresponds to the replacement of 
supertrace by ordinary trace. Our conventions for supertrace (and superdeterminant) are defined 



in Appendix IX 
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In the proof of the Polyakov-Wiegmann identity we exploit the cyclic property of the trace. For 



supermatrices, it is the supertrace which has this property (see Appendix IX). The form of the 
Polyakov-Wiegmann identity for supermatrices is the same as that for ordinary matrices, providing 
we replace trace by supertrace. 



16 



The action ( |5.7|) is invariant under the semi-local transformation 

g(z,z)^Q(z)g(z,z)n- 1 (z) (5.8) 

where Q(z) and Q(z) are arbitrary elements of G. This invariance is made manifest by the 
Polyakov-Wiegmann identity ( |5.4|) with A 2 = 47r/fc, 

W[ab] = W[a] + W[b] + — [ d 2 £ str (a -1 da bdb' 1 ). (5.9) 

2n J 

One may consider promoting this semi-local symmetry to a true local symmetry, with the 
introduction of auxiliary gauge fields. In particular, the action 

W[g J h } h] = W[h- 1 gh]-W[h- 1 h] (5.10) 

is clearly invariant under the combined (gauge) transformations 

g — > X(z, z)(yfA _1 (z, z), h—>\(z,z)h, h^X(z,z)h. (5-11) 

Applying the Polyakov-Wiegmann identity ( |5.9| ) and defining the following gauge fields 

A = hdh-\ A = hdh- 1 (5.12) 

one may rewrite the action ( |5.10| ) in the following form, 



W[g,A,A}=W[g} + ^ j d 2 fstr (Ag^dg - Adgg~ x + Ag^Ag - AA) (5.13) 
= W\9\ + rf ^ str ( aJ a + AJ A ) (5.14) 

where we have introduced the gauge invariant generalizations of the conserved currents ( |5.6D : 

J A = g-\B + A)g, J A = -(d + A)gg- 1 . (5.15) 

The gauge fields A and A are non-propagating, and play the role of Lagrange multipliers 
forcing certain currents to vanish. Choosing the gauge fields so that the currents associated 
with the group H are set to zero, we may describe the coset space G/H. The action (|5.13|) 
will be our staring point in the next section. 



VI. DECOUPLING AT THE KAC MOODY POINT 

In order to see how the SL(2, C)/SU(2) WZNW may emerge as an independent subsector 
of the Zirnbauer model (at the Kac-Moody point) or some other 'grand' supermanifold 
theory, we parametrize the supergroup manifold using the Gauss decomposition. The use of 
such decompositions is well established in the the free-field approach to WZNW models (see, 
for example |39| and references therein). In particular, such an approach has been followed 
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in the study of the gl(n|n) current algebras and associated topological field theories flO . 
We consider the gauged WZNW modeQ, 

W[g, A, A] = W[g] + A J str (Ag^dg - Adgg~ l + Ag~ l Ag - AA) (6.1) 

in which the supergroup element, g G SL(2|2), admits the following Gauss decomposition 
(see equation |A7|) 

9 = 'r, i=(*°)(°?)(U) (6-2) 



x A y \0 b) \o u / 

where $ G C, A and x are arbitrary 2x2 Grassmann-odd matrices, and a and b are arbitrary 
2x2 unimodular matrices (det(a) = det(fe) = 1) with Grassmann-even entries. One may 
recover the Zirnbauer's 'base manifold' by taking a G SL(2, C)/SU(2) and b G SU(2). In 
order to set the gl(l) currents equal to zero, and thus describe the PSL(2|2) coset, we choose 
(see equation | A 1 0|) 



Applying the Polyakov-Wiegmann identity (|5.9| ) to the first term of (16.11) , 

W[g] = W[e*\ + W[y\ + — f 9$ str (-fd^ 1 ) = W[y] (6.4) 

where we have used the fact that W [e*] = and that str (jdj" 1 ) = 0. Repeated application 
of the Polyakov-Wiegmann identity ( |5.9| ) shows that 

W[j) = W[a] - W[b] + — [ tr(b- l d\adx)- (6.5) 

2n J 



Using the fact that str(E • ■ ■ ) = tr(- • ■ ) for arbitrary arguments, the second term in ( pTT 
may be written 



k 

2^ 



Straightforward matrix manipulation of 7, together with the fact that tr(a 1 da) = 0, 
tr(<9aa -1 ) = 0, and likewise for b, yields 

trfr -1 ^) = 2tr (b^dXax) , (6.7) 

tr(d77~ x ) = 2tr (b~ l \adx) , (6.8) 

tr(7 -1 E7) = -4tr (b^Xax) ■ (6.9) 
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Throughout this section f = J d 2 ^. 
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Combining the above results we may write 
W[g,A,A} = W[a\-W[b] + — [ [tr (b' 1 {8 - ft)\a(d + /i)y) + 2/i<9$ - 2/29$l . (6.10) 

27T J 

It is convenient to introduce new fermionic fields A' and x' via 

A = 6V, x = a~V- (6-11) 

Since a and 6 are unimodular matrices, the corresponding Jacobian is unity. The trace term 
becomes 

tr [(8 - jl + b~ x 8b) \'(d + fi + ada- 1 ) X '] = tr [(8 + A b )X'{d + A a ) X '] (6.12) 

in which we have reparametrized, /i = da, and p, = 8(3, and introduced the gauge fields 

A b = (e /3 r 1 )9(e- /3 6), A a = {e' a a)d{e a a - 1 ). (6.13) 

We make a further change of fermionic variables 

\"={B + A b )\>, X "=(d + A a )x' (6.14) 

and note that one may express the corresponding Jacobian in terms of the WZNW action 
on the matrices a and b (not b^ 1 ) respectively [38, 41|. Recalling that the fermionic fields 
are 2x2 matrices, this leads to a level shift of 2 in the WZNW models defined over a and b. 
Such level shifts were also encountered in [ 34^ 40|1 . One obtains the following action for the 
PSL(2|2) WZNW model, 

M[PSL(2|2)] = (2 + k)l[a] + (2 - k)l[b] + ■■■ (6.15) 

in which the ellipsis indicates terms independent of a and b; these include the contributions 
of the bosonic fields a and f3, the fermions A" and x" an d ghosts. In the case of the Zirnbauer 
'base manifold', a gSL(2, C)/SU(2) and b eSU(2). The level shifts are important here, and 
the resulting contribution to the central charge from these (renormalized) WZNW models 
is level- independent: 

3(fc + 2) 3(fc-2) 
C= (fc + 2)-2 + (A:-2) + 2 =6 ((U6) 

In particular, one observes that the SL(2, C)/SU(2) WZNW model emerges as an indepen- 
dent subsector of the theory. 



VII. CORRELATION FUNCTIONS OF LOCAL DENSITY OF STATES 

We have provided evidence to suggest that the correlation functions of the normalized 
local DOS may be obtained from the critical SL(2, C)/SU(2) WZNW model, in which the 
level k is fixed by the known scaling dimensions (|2.11| ). This model has been studied before, 



even in the context of the theory of disorder . From the preceding sections, we have seen 
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how this theory may be embedded, for example, into a larger theory such as the PSL(2|2) 
model (at the Kac-Moody point). Taking into account the level shifts induced in this 
embedding one may write the action for the SL(2, C)/SU(2) WZNW model ( 4.15 ) in the 
following form 

S = <*±9 J d 2 x [48686 + e^dfidfi*] (7.1) 



A comparison between the decomposition of SL(2, C)/SU(2) appearing in ( [4. 14|) and the 
Gauss decomposition of SL(2,R) 



7 I J \0 e-"J \0 1 

where 7, 6, <fi G R is convenient at this stage. In particular, one may obtain the parametriza- 
tion of SL(2, C)/SU(2) from that of SL(2,R), by first complexifying the real parameters 
7 and ip, and then making the identification 7* = ip. In this way, one may obtain the 
SL(2, C) /SU(2) correlation functions from the (much studied) SL(2, R) WZNW modeQ one 
treats 7 and ip as independent real fields and performs the appropriate continuation wherever 
necessary. 

The action ( |7. 1| ) , is characterized by a single coupling constant k, which determines the 
scaling dimensions: 

d 3 = -2j{j + l)/k (7.3) 

There are no restrictions on the value of k, and the values of j are related to representations 
of the SL(2,R) group (see Appendix |A~^). 

As we have discussed in Section [TTI| , the available information about the local DOS 
appears to be self-contradictory. On the one hand, the fusion rules for [^(r)] 9 suggest 
that these are operators belonging to the sl(2,R) representation with j = —q and with 
scaling dimensions given by (|7.3| ). On the other hand, the two-point correlation function in 
the strip geometry is given by an integral over continuous representations of sl(2,R), with 



j = — 1/2 + ip/2. As we have demonstrated in Section |T|, this paradox is resolved if, in 
fact, the two-point function of the local DOS is given as the limit of a four-point function 
in our field theory: 

0,(1,2) = lim <M(l)M(2)*o(3)*o(4)> (7.4) 

|r3|=0,|r4|— >-oo 

The four-point correlation function appearing in this relation (which we shall subse- 
quently denote by G q (l, 2, 3, 4)) is determined by the Knizhnik-Zamolodchikov equations 



15 Despite the fact that the SL(2, R) model cannot be defined as an Euclidean path integral, it 
can be defined algebraically. Following this procedure one can derive the Knizhnik- Zamolod- 
chikov equation for the correlation functions. This equation is not different from the corresponding 
equation for the SL(2, C)/SU(2) model. 
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for the SL(2,R) WZNW model; these equations have been derived in |37|,^2|. The operators 
in (|7.4j) are matrix elements of infinite-dimensional matrices realizing particular representa- 
tions of the group (information about representation theory of sl(2,R) is given in Appendix 
4). Instead of working with these particular matrix elements, it is more convenient to calcu- 
late the correlation function between all possible tensors belonging to representations with 
j = 0, j = —q. We emphasise that one does not need to worry about exact definition of 
\l/o because the Knizhnik-Zamolodchikov equations together with the asymptotics of the 
solution automatically select the necessary primary fields. By that token standing in 
the correlation function is the very primary field with zero scaling dimension we need. It 
is necessary to mention here that there may be several fields with zero dimensions with 
different symmetry properties and unusual OPE ^3J leading to other solutions of Knizhnik- 
Zamolodchikov equations for correlation function ( |7.4| ). However these solutions are excluded 
since they do not lead to the desirable quasi-one-dimensional asymptotics. 

We introduce the fields p q (z,y), ^/ (z,y) as described in Appendix 4. Now each point in 
( [7.4| ) is characterized by two complex coordinates z and y. Invariance under the projective 
transformations of the plane[^], together with the sl(2,R) Ward identities, restricts the four- 
point function to have the form, 

g,(l,2,3,4)= 1 z = ^ t= y^L } (7 . 5) 

\Zl2r hq \yi2r q Z 31 Z i2 2/312/42 

where we have used the fact that [p q ] has conformal dimension h q (= d q /2), and sl(2,R) spin 
j = —q, whereas \&o has zero conformal dimension and j = 0. By definition, the function 
is invariant under the permutation of Z\ and z 2 which implies 



F q (z,z;t,t)=F q (l/z,l/z,l/t,l/t). (7.6) 
In the limit, — > 0, corresponding to the fusion of the local DOS, one obtains 
([^](1)[/](2)* Q (3)* Q (4)> - \z 12 \~ Ah ^ ([^](2)* Q (3)* (4)> 

-yOO 
2q 



where C 2q is a structure constant of the conformal field theory. This fixes the asymptotics 
of T q {z) at 1 — z, 1 — y — > 0: 

F q (z,y)~(l-z) h *>y 2 * (7.8) 

We may explore different asymptotics of J- q {z) in the strip geometry. Under the transfor- 
mation z = exp(w/R) the point 3 goes to — oo and the point 4 goes to +oo. Considering 
only the holomorphic part and omitting the y-dependence, one obtains 

£,(1,2,3,4) = [2Rsinh(w 12 /2R)]- 2h «F q [exp(w 12 /R)\ (7.9) 
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Namely transformations of the form w = (az + b)/(cz + d), where ad — be = 1. 
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Thus in the strip geometry the correlation function is translationally invariant. The limit 
Rewi2 3> R is governed by the asymptotics of T q {z) as z — > oo and the limit Rew 2 i R 
is governed by the asymptotics as z — > 0. Both limits are related by the crossing invariance 
condition (|7.6|) . 

Thus, the fusion rule ( |2.9| ) and the expansion in continuous representations ( 2.12|) appear 
in different channels of the four-point function, thus resolving the apparent paradox. 

According to |J7|,^2| the Knizhnik-Zamolodchikov equation for model ( |7. 1| ) is given by 



(v«) a 



o 2 



dyidy. 



+2( ^lr Ji 4)~ 2 ^ 



G(l---4) = (7.10) 



as may be seen by combining equations ( |A30| ) and ( |A45| ) together with a suitable redefinition 
of k. Substituting (|7.5| ) into this equation, and using the global conformal invariance to map 
(zi, z 2 , z 3 , Za) — > (1, z, 0, oo) (and similarly for y using the sl(2,R) Ward identities) one arrives 
at the following partial differential equation for T{z, t) 



p (t-z)t(t-l) + 
z{z-l) 



(1-t) 2 t(2q-t) 



z-1 



+ kd z T = 



(7.11) 



iFrom equation (|7.8|) we conclude that in the 'discrete' channel \z — 1| <C 1, T — (z 
^-) h2q f{t)- Equation ( [7.11]) may now be written as 



(l-t)%[td t f] + kh 2q f = 0. 

The solution of this equation takes the form 

f[t = 1/(1 -x)] = F(2q,l-2q,l;x), 

= [d x ] 2 "- 1 [^(l-a:)- 1 ], 

2q-l 

= (2q - - 1) 2 « + J2(i- t- i y k k\c k £i 

k=0 



(7.12) 



(7.13) 



Since equation (|7.12|) is invariant under the transformation t — > 1/t the second linearly 
independent solution is 

f(t) = /(1/t) (7.14) 
and the entire crossing invariant solution at \z — 1\ « 1 is given by 
F(z,z;t,t) = \1 - z\ 2h ^{A[f(t)f(t) + f (1/t) f (1/t)} + B[f(t)f(l/t) + f(t)f(l/t)\} (7-15) 



Following equation ( |3.6| ) at \z\ <C 1 (the 'continuous' channel) we look for the solution 
in the form of a linear combination of power law solutions: 



T = J df,(p)\z\ 2D ^ k F p (t,t) 



(7.16) 
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where 

D(p) = (1 + p 2 )/A - kh q = {q- 1/2) 2 + p 2 /A. (7.17) 
The function F p obeys the equation: 

t 2 (l - t)d 2 t F + t(2q - t)d t F + D(p)F = 0. (7.18) 
After the change of variables t — 1/r one obtains the hypergeometric equation: 

r(l - r)d 2 F + [1 - (2 - 2q)r]d t F - D(p)F = (7.19) 
There are two independent solutions of this equation: 

FW(t) = F(j,j*,V,l/t), 

F^(t) = \ntF( Jj3 \l;l/t) + R(l/t), (7.20) 

where j = 1/2 — q + ip/2 and R(x) is regular as x — > 0. Thus, the general solution ( [7.1 6| ) at 
small |z| is 

T= I d^p)\z\ 2D ^l k C ab {p)Fj?\t)F^\t) (7.21) 



where the matrix C a b (a>,b = 1,2) is determined by the crossing symmetry ( |3.3j ). At the 
moment we do not have a complete solution of equation ( [7.1 1|) which prevents us from 
determining all the constants in our equations and OPEs. It appears, however, that the 
solution is possible. Here we shall just outline the method leaving the detailed analysis for 
our next publication. 

The general solution contains the parameters i/i which replace matrix indices for operators 
belonging to infinite dimensional representations. Presumably the local DOS corresponds to 
certain matrix elements. In order to identify p q it is instructive to take a closer look at the 
solution on the strip. Combining equation (|7.9| ) with equation ( [7.21|) we obtain the following 



expression for the asymptotics of the two-point correlation function of operators p q (w,y): 

p*(w uyi )p«(w 2 ,y 2 ) = \y 12 \-* q J dfi(p) exp [-(1 + p 2 )\r 12 \/2kR] C ab (p)F a (t)F b (i) (7.22) 

where w = r+ia (0 < a < 2ttR). In the limit |r 12 | ^> R the prefactor of the integral becomes 
unity. One recovers the correlation function ( p.!5[ ) leaving in the expansion of (|7.22| ) in t, i 
only the term t~ 2q i 1 ~ 2q . 

We use this fact to determine which matrix element represents p q . Therefore let us 
assume that this is a matrix element of p q (y) containing y n . In order to extract it from 
the y-dependent correlation function one has to integrate it with the weight y~ x ~ n over a 
contour in the y-plane surrounding zero. Then we have 

/ f[dy i y^- n y 2 1 - n y^ q Ht) = 

J i=l 

d^dU^r^klt^ J dtr 2n - 2q+2 Ht) j d * {l _Jr x _ t) 2 ~ 

J drF(j.f,l;r)-^ 



-r) 



T r 2(q+n-l) 
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where we have performed the transformation 

2/32 Z/il/32 - ty 2 yai 
x = — , 1/4 = . 

2/31 2/32 - tffyi 

Notice that we integrate over y%, y& with unit weight. In order to recover the term r 2<2_1 in the 
expansion of the hypergeometric function corresponding to the necessary matrix element, 
one has to put n = —2q. Thus we identify p q with the highest weight of the j = —2q 
representation. 

The Knizhnik-Zamolodchikov equation (|7.11| ) can be solved by the method suggested 
in [ 4~2| . Let us make the change of variables: 

T{z, t) = t a (t - l)\t - z) c z\z - l) e Y(t, z) (7.23) 



where 



a = -(l + /e)/2, b = q - (k + l)/2, c = 1/2 - q, 

d = -(2q-l)/2k, e = (2q-l) 2 /2k. (7.24) 



In terms of Y, equation ( |7.11| ) becomes 



+ [h^t - z)~ 2 + h 2 r 2 + h 3 (t - I)' 2 - n[t(t - l)]- 1 ]^ (7.25) 



where 

K = h 1 + h 2 + h 3 -h 4 -(3k + 2)/4 (7.26) 

This equation coincides with the equation for the 5-point function of vertex operators in the 
Liouville theory 



St= d x 



k_ 

4tt 



{d^f + 2{k + 1)1Z4> + r/exp(-20) 



(7.27) 



(1Z is the Riemann curvature of the world sheet) with the central charge equal to 

C = 1 + 6(Vk + l/Vkf (7.28) 
The operators of the WZNW model are identified as vertex operators of the Liouville theory: 

V n = exp(n0) (7.29) 

with certain ns. 

This establishes a link between the Liouville theory and the theory of the local DOS at 
the plateau transitions in the IQHE. 
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VIII. CONCLUSIONS 



We summarize here the results of this paper: 

1. We have described a general procedure for calculating correlation functions of local 
DOS at the plateau transitions in the IQHE. This procedure requires insertions of 
additional vacuum operators into any correlation function, so as to modify the ground 
state. 

2. We have provided arguments, that the correlation functions of the local DOS are 
appropriately described in terms of the SL(2, C)/SU(2) WZNW model at the usual 
Kac-Moody point and with a suitably chosen level k. The comparison with the nu- 
merics gives 6 < k < 8. In this model we have identified the operators corresponding 
to the properly normalized local DOS, and calculated their correlation functions. 

3. We have demonstrated that the SL(2, C)/SU(2) WZNW model may emerge as an in- 
dependent subsector of some supergroup theory, and may indeed describe a disordered 
system. 

There are many unresolved problems which require further study. A number of these are 
outlined below. 

The model contains a free parameter k which we fix by comparison with the known 
scaling dimensions. It is not clear at the moment what mechanism is responsible for pinning 
k to this particular value. The example given in Appendix |A6| demonstrates the kind of 



surprises one can expect in dealing with theories on non-compact manifolds. 

We have studied in some detail only the correlation functions and scaling dimensions 
of the local DOS. We are not in a position to discuss the correlation function exponent v 
which is, presumably, determined by an operator which does not belong to the non-compact 
bosonic sector of the theory. 

Even in the non-compact bosonic sector of the theory there are important unsolved prob- 
lems. For example, we have not yet solved the Knizhnik-Zamolodchikov equation (|7. 1 1|) 
which determines the two-point correlation function of the local DOS. In subsequent publi- 
cations we shall describe a procedure for calculation of the multi-point correlation functions. 
This will allow us to explain the termination of the multifractal spectrum. As we have 
mentioned in Section |11 C| the multifractal spectrum of the local DOS terminates at q ~ 2.5. 
If the SL(2,C)/SU(2) WZNW model is the correct model for the local DOS, it must possess 



some termination mechanism. This may be a mechanism similar to the one described in |3T 
in the context of the theory of Dirac fermions with gauge potential disorder. The latter 
mechanism was related to the logarithmic nature of the theory. 
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APPENDIX: 
1. Supergroups 

In this section we provide a brief introduction to supergroups. For the reader interested 



in further details we refer them to the book |12[ and references therein. For the purposes of 
this paper we consider matrices of the form 

M = (" f) (Al) 



P b 

where a and h are respectively n x n and m x m matrices containing even elements of a 
Grassmann algebra, and a and p are respectively n x m and m x n matrices consisting of 
odd elements of a Grassmann algebra. Such matrices are called supermatrices. The set of 
such complex (respectively real) square supermatrices is denoted by M(m|n; C) (respectively 
M(m\n; M)). The supertrace of M is defined as 

str(M) = tr(a) - tr(6) (A2) 

where the symbol tr stands for the conventional trace. This definition provides the invariance 
under cyclic permutations: 

str(MiM 2 • ■ • M n ) = str(M n Mi • ■ • Af n _i) (A3) 

for arbitrary supermatrices Mi,-- - ,M n . The superdeterminant (or Berezinian) of M is 
defined as 

, /„,s det(a — ab' 1 p) , . i% 

SdCt(M) = detW (A4) 

This definition provides the factorization property of the superdeterminant: 

sdet(M!M 2 • --Mn) = sdet(M!)sdet(M 2 ) • • • sdet(M n ). (A5) 

The general linear supergroup GL(m|n; C) (respectively GL(m|n;R)) is the supergroup 
of invertible complex (respectively real) supermatrices, the group composition law being the 
product of supermatrices. In particular, an arbitrary element g G GL(n\n) admits the Gauss 
decomposition [40 



11 0\ (A 0\ ft x 



(a nJ(o b){o i) < A6 ' 

where A and B (A and x) are arbitrary nxn Grassmann-even (respectively Grassmann-odd) 
matrices. 
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The special linear supergroup SL(m|n;C) is the subsupergroup of supermatrices M G 
GL(m|n; C) such that sdet(M) = 1. This requirement forces det(A) = det(B) in the Gauss 
decomposition (|A6|). Thus, an arbitrary element g G SL(n\n; C), may be decomposed as 

— <*> (a S) (o t) (S J) (AT) 

where $ G C, and a and 6 are arbitrary n x n unimodular matrices (det(a) = det(6) = 1) 
with Grassmann-even entries. 



2. Superalgebras 

In addition to supergroups, one may also introduce the notion of a Lie superalgebra, in 
which the Lie bracket is replaced by a generalised bracket (commutator/anticommutator) 
that depends on whether the generators considered are 'bosonic' or 'fermionic'. A particu- 



larly useful dictionary of Lie superalgebras has been compiled in |44[ . 

The superalgebra gl(2|2) is generated by sixteen independent matrices, which in a suitable 
basis may be chosen as the twelve off-diagonal matrices consisting of a single entry of unity, 
together with the matrices 



13 



where a 3 = diag(l, — 1), and II denotes the 2x2 unit matrix. The Lie algebra sl(2|2) is 
generated by supertraceless matrices, and may be obtained from the above representation 
by removing X. An arbitrary sl(2|2) valued current may be expanded as 



ii 



J{z) =i{z)I + a{z)E + ^2t i {z)T i . (A9) 

i=i 

The Lie algebra sl(2|2) contains a non-trivial centre; the unit matrix / commutes with all the 
generators. The algebra psl(2|2) is obtained by removing this generator. One may isolate 
the contribution to the sl(2|2) current arising from the identity component by the following 
relation 

i(z) = ^str(SJ) = -trJ (A10) 



3. Scaling Dimensions 

We take a <i-dimensional sample of material of side L and divide it into = (L/l) d boxes 
of side /. The probability of finding an electron in box i (the so called box-probability) is 
given by 




(All) 
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The average value of the box probability scales in the following manner 

N / i \ d 



8=1 V 7 



(A12) 



for all normalized wavefunctions (^2 i Pi = 1) and is not useful in distinguishing between 
localized, extended, and critical wavefunctions. We are led to consider the scaling of the 
moments of the box-probabilities 



(P q ) 



which serves to define new exponents r q . Consistency with ( |A12j ) requires^ 

7~i = To = —d. 



(A13) 



(A14) 



It is often advantageous (especially in numerical simulations) to introduce the Legendre 
transform, /, of r(q), 



f(a g ) = qa q - r ( 



where 



It follows that 



a q = 



dr q 
dq 



df(a) 
da 



(A15) 



(A16) 



(A17) 



Thus, the function f(a) has a maximum at q — 0. Using ( |A14| ) and ( A15|) , 

/(a) max = f(ao) = d. (A18) 
The slope of f(a) is unity at the value of a corresponding to q — 1, 

df(a) 



da 



1. 



9=1 



Combining (|A14 ) and ( A15 ) we also find that, 

f{a l ) = a 1 



(A19) 



(A20) 



17 For a uniform electron distribution (-P 9 ) uniform = j? YliLi i~k) 9 ~ (t)^ anc ^ one ^ n ds that 
T q = {q — \)d. It is conventional to define generalised dimensions, D q , such that r q = (q — l)D q . 

The scaling relation ( A13| ) may equivalents be written (P q ) ~ {{) d+{q l)Dq ■ We note that ( [All) 
implies D = d. 
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The constraints ( A18|) , ( |A19|) and ( |A20|) allow us to write a parabolic approximation in terms 
of one parameter, o (the position of the maximum)^] 



f(a) = d 



(a - g ) 2 
4(o — d) 



(A21) 



This parabolic approximation together with the definitions ( |A15|) and (|A16|) gives rise to 
the relation 



Tq - a 0) 



dr 

q— r = d — 

dq 4(ao — d) 

This equation is solved exactly by the polynomial 

r q = -d + a q - («o - d)q 2 . 



Recalling the scaling relation (|A13| ) one obtains 



(P q ) 



L 



aoq-(ao-d)q 2 



The local density of states is given by 



p = 5~ l P ~ L d P. 



Hence 



(A22) 



(A23) 



(A24) 



(A25) 



(P q ) 



L{a -d)q(l-q) ' 

Numerical simulations for d = 2 give «o ~ 2.28. 



(A26) 



4. Representations of the SL(2,R) Group 

The Lie algebra sl(2) is generated by three independent traceless matrices, which, in the 
spin basis may be taken as 

■r-02). '-Ki-O. "-(Si)- (A27) 

Their commutation relations read 

[J + , J~] = 2J°, [J°, J±] = ±J ± J (A28) 



18 Assume f(a) = d — A(a — ao) 2 where A is a constant to be determined; this takes into account 
(HKWj ). Condition ( gig ) implies ai = ao — 1/(2A). Substituting this into the equation following 
from ([A20p , namely a\ = d — A(a\ — ao) 2 , determines A. 
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and should be familiar from the theory of angular momentum. It is readily seen that the 
differential operators 

J- = ^~, J° = y^--j> J + = 2jy-y 2 ^- (A29) 
ay ay ay 



obey the same commutation relations ( |A28| ) when acting on the space of differentiable func- 
tions. 

In this differential realisation, the monomial yi +m plays the role of the state vector \ j, m) . 
One observes that d/dy acts as a kind of lowering operator, and y acts as a kind of raising 
operator. It is convenient to modify our basis slightly by defining t +1 = —J + , t° = J°, 
t^ 1 = J~. This enables us to write the generators in the following compact form 

t' = / +1) ^-(/ + i)jV; / = -i,o,+i (A30) 

Their commutation relations read 

[t + ,r] = -2t°, [t ,**] = ±t ± (A31) 
That is to say we have structure constants, 

/° + + = 1 = -/°-_, / + "o = -2 (A32) 
The Killing form reads (with the ordering +, 0, — , and c v = 2) 

/0 2\ /0 §\ 

rj ab = -1 , Vah = -1 (A33) 
\2 0/ V| 0/ 

More generally, one can construct representations defining their action on functions f(y): 

Af(y) = \cy + d\ 2 i sign (cy + d) 2 J (^^) (A34) 

\cy + dj 

where ad — cb = 1. 

The operator p q discussed in the text belongs to the j = —q representation and is 
annihilated by t~\ 

t+ = _ q p q = 0, P H- = (A35) 
In other words one may perform the following expansion 

P q (y,y)= £ v-^V t *£? ) (A36) 

k,k=0 

In this representation t 3 is diagonal : 

t 3 $ krk = -(q + k)$ krk (A37) 
The operator $ belongs to the representation where t 3 cannot be diagonalized: 

t + $o = ti$ , $ r = (i$ (A38) 
For more information see, for example 
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5. The Knizhnik-Zamolodchikov Equation 

The stress-energy tensor of WZNW models is proportional to the scalar product of 
currents. Therefore the Virasoro generator with n = — 1 acting on any primary field gives 



k + c v 

We consider the insertion of the zero vector 



-2 

k + c 



\x) 



L-i + 



k + 



= o 



(A39) 



(A40) 



inside the correlation function of a set of primary fields. We note that the insertion of the 
operator J" in the correlator can be expressed as 



yizx) ■ ■ ■ {J°Li(j)i){Zi) ■ ■ ■ (j) n {z n )) 



dz 



2ni /„. z — Zi 



(J a {z)(f> 1 (Zl)---<f>n^n)) 



(A41) 



r\ 7 ^ 

^ ^-^-V ^ '-^ (A42) 



E 



HOO • • -(j)n(z n )) 



(A43) 



Therefore 

(01 («i) • ~x(?i) ■ ■■ ( Pn{z n )) 



k + c. 



E 



(A44) 



and by construction this must vanish: 



A. - 



E 



L (Zl) ■ • -4>n(z n )) = 0. 



(A45) 



This is the Knizhnik-Zamolodchikov equation. The solutions to this equation are the corre- 
lation functions of primary fields. 



6. BRST Invariance and A = Conformally Invariant Deformation 

Let us consider a WZNW model at level — k^\ deformed by its kinetic term: 



8(9) — SwZNwid 



-1\ 



d 2 zVl(z, z) 



(A46) 



19 The notation is choosen in such a way that for k > the action for the SL(iV, C)/SU(A r ) model 
is positive. 
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2iri( 


— z 


d( 


1 


2m C 


— z 



where 

Q(z, z) = J a {z)J\z)(t) ab {z, z), (A47) 

and 

J = hg' l dg, J = l -kdgg-\ ^ = tr^V^ 6 . (A48) 

The perturbation operator fl possesses a number of interesting properties. Let us com- 
pute the commutator of the Kac-Moody current J a (z) with the operator Q. Denoting its 
chiral part O(z) we get 

[r(y),o(z)]= /^^{[J a (y),JWW + AC)[^(y),0 b W]} 

W b J c (0<j> b (z)5(y, C) + l -k<t> a {z)5'{y, Q + /fJ b (C)0 c W%, *)} 

fab £cb i fab tbc 

!f^L<p d (z)5(y, + l -knz)5\y, Q + ^<t> d {z)5{y, z) 
Q- z I Q- z 

+ <fi ^T— [f^ cd (z)6(y,() + f^ cb {z)5{y,z)}. (A49) 
Here 

m c \ z ) = : J c {z)(t) b {z) : . (A50) 
By taking contour integrals, we obtain 

[J\y), O(z)] = {- l -k + c v )<f> a (z)5'(y, z). (A51) 
Thus for the entire operator f2 we get find 

\J a (z),n(w,w)] = (-^k + c v ^j J b (w)<p ab (w,w)5'(z,w). (A52) 

For the SL(iV, C)/SU(A r ) model where cy = N this commutator vanishes when k = 2N. 
(as in the main text the level k is defined such that the central charge is given by C = 
k(N 2 — l)/(k — N)). 

In other words, when k = 2N, the operator Q is invariant under the Kac-Moody symme- 
try. Since the Virasoro operators are quadratic combinations of the Kac-Moody currents, the 
operator Q is automatically invariant under the conformal transformations. In particular, 

[L ,n] = Afi = 0. (A53) 

Indeed, 

A(k = 2N) = 0. (A54) 
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The operator £7 forms a closed OPE algebra 



n(z)£l(0) = fi(0). (A55) 

The constant term on the right hand side is absent because the operator Q has zero norm; 
this follows from equation ( |A52| ). 

Thus, a the SL(JV, C)/SU(JV) WZNW model at level k = 2N perturbed by the operator 
fl must be conformally invariant for arbitrary parameter e. This can be seen in the following 
way: The property (|A55|) implies that away from the conformal point the following relation 



holds 

dT = d&, (A56) 

where T and are the components of the stress-energy tensor. Moreover, the trace is given 
as follows 

e = fa n, (A57) 

where fa is the renormalization group beta-function of the coupling e. Since Q is not a 
marginal operator, fa ^ 0. However, the following is true 

06= [L_!,e] =#)[£-!, ft]- (A58) 

Since [L n , Q] — for any n, we arrive at 

dT = 0, (A59) 

i.e. the stress-energy tensor component T is still a holomorphic function even away from 
the critical point. This means that the perturbed theory remains conformal for an arbitrary 
perturbation. This also means that the most stable point in this case is not k = 1 as for the 
SU(N) WZNW model, but k = 2N. 

The described possibility is not realized for the PSL(2|2) model where cy = 0. Though 
the scaling dimension of the Q operator in this case is equal to 2 and it may appear being 
exactly marginal, the fact that it does not commute with J's may be an indication that 
the line of critical points obtained by studying perturbative series @, [35| is destroyed by 
non-perturbative effects. This is an open interesting question. 
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